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ABSTRACT 

We present an analysis of the globular cluster NGC 2419, using a polytropic model in Modified 
Newtonian Dynamics (MOND) to reproduce recently published high quality data of the structure and 
kinematics of the system. We show that a specific MOND polytropic model of NGC 2419 suggested 
by a previous study can be completely ruled out by the data. Furthermore, the highest likelihood fit 
polytrope in MOND is a substantially worse model (by a factor of ~ 5000) than a Newtonian Michie 
model studied in an earlier contribution by Ibata ct al. (2011). We conclude that the structure and 
dynamics of NGC 2419 favor Newtonian dynamics and do indeed challenge the MOND theory. 
Subject headings: dark matter - galaxies: clusters: individual (NGC 2419) - gravitation 



1. INTRODUCTION 

The well-known discrepancy between the luminous 
components of galaxies and their dynamics has been in- 
terpreted as evidence that vast quantities of some un- 
known type of dark (non-luminous) matter surround 
these celestial structures. Although this dark matter has 
not yet been directly detected, its existence is consistent 
with other astrophysical analyses, including the dynam- 
ics of galaxy cluste rs and the formation of large-scale 
structure (see, e.g., Springel et al. 2006). However, a 
ver y interest i ng alt ernative to this possibility, proposed 



by Milgrom| (1983), is that our theory of gravity is in- 
correct or incomplete. According to this idea, General 
Relativity (or the Newtonian approximation to that the- 
ory) in the low-acceleration regime below a characteristic 
value ao(~ 1.2 x 10~ 8 cms~ 2 ) under-predicts the actual 
acceleration. 

This Modified Newtonian Dynamics (MOND) theory 
has succeeded in passing numerous observational tests 
over the almost three decades since it was first proposed 
(for instance, it is able to fi t the rotation curves of low 
surface brightness galaxies, McGaugh & dc Blok |199"8l 
and tidal dwarf galaxie s, IGentile et al.||2007 [ and can 
fit gravitational lenses, Shan et al. 2008p . Although 
there are outstanding problems with the predictions of 
the M OND theory (the growth o f cosmological perturba- 
tions, [DodeIson~^JLiguori 2006, the offset b etween lens- 
ing m ass and baryons in t he Bullet Cluster, |Clowe et al. 
20061 Sol ar System tests, |Milgrom|[l983l |Sereno fc JeF 
zer||2006 dynamical friction in dwarf galaxie s^ |Ciotti & 
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to list a few) , the more widely-accepted Cold Dark Mat- 
ter (CDM) paradigm has deep-set issues of its own (the 
missing satellite galaxies, the possible non-existence of 
dark matter c usps and the high angular m omentum of 
galactic disks, |Binney]|2004l [Primack |2009[ p] Determin- 
ing which of the two competing theories of gravity repre- 
sents reality remains a very important and fundamental 
task. 

Recent improvements in the precision and multiplex- 
ing capabilities of spectrographs have led several teams 
to study large samples of stars in Galactic globular clus- 
ters as a me ans to test MOND on the scale of parsecs up 



to ~ 100 pc (Scarpa et al.|2003| |Haghi et al.||2009| Jordi 



eF^[2009| [Gentile et al.||2010| |Lane et ai.||20TU| f. lna 
recent contribution to this effort, our team examined the 
distant halo globular cluster NGC 2419, re-measuring the 
surface brightness profile derived from HST/ACS, Sub- 
aru/Suprimecam and CFHT/MegaCam images, and ob- 
taining high accuracy radial veloc ity measurements with 
the Keck/DEIMOS spectrograph ( jlbataet al.|2011[ here- 
after Paper I). We argued that this globular cluster is by 
far the best globular cluster target to test MOND with 
current instrumentation. Most importantly, it lies far 
enough from the Milky Way so as to relatively isolated, 
and so suffers little external acceleration. Furthermore, 
it is massive, so it offers a reasonable number of targets 
within reach of available high-resolution spectrographs, 
its velocity dispersion is large enough to be easily mea- 
surable, and it is spatially extended, meaning that a sub- 
stantial fraction of the cluster lies in the low-acceleration 
regime where MOND differs from Newtonian gravity. 

Outside of the very central regions, the cluster is 
circularly-symmetric in projection, furthermore, the clus- 
ter does not show evidence of significant rotation (Pa- 
per I). These facts justify using spherical models to fit 
the system. The analysis presented in Paper I first ex- 
amined Michie models, which are an extension to King 
models that allow for anisotropy in the orbits of the con- 
stituent stars of a spherical cluster. We found that a 
Michie model in Newtonian gravity (with a core radius 

6 Plausible simulated solutions for these problems with CDM 
have only recently been presented, yet they remain untested obser- 
vationally. 
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of r c = 11.7pc, a tidal radius of r t = 331 pc, a central 
velocity dispersion parameter of <r — 6.36 kms -1 , and 
an anisotropy radius of r a = 1.5 x r c ) gives an excel- 
lent representation of the data, whereas Michie models 
in MOND have significantly lower likelihood (the best 
model being a factor of more than 40000 worse than the 
best model in Newtonian gravity). 

We then expanded our analysis to examine a wider 
range of models, using a Markov-Chain Monte Carlo 
(MCMC) scheme to generate kinematic models consis- 
tent with the Jeans equation. The best such model in 
MOND was again found to be strongly disfavored com- 
pared to the best Michie model in Newtonian gravity. 

However, these con clusions from Pa per I were recently 



strongly criticized by Sanders ( 20lT] hereafter Sll), on 
the grounds that other models could reproduce the sys- 
tem well in MOND. The MCMC procedure mentioned 
above that we implemented in Paper I examined one mil- 
lion different kinematic models. Nevertheless, it is cer- 
tainly possible that a good model in MOND was missed 
by the MCMC routine. Indeed, Sll claims to have found 
a model that provides a counter-example refuting our 
conclusions, proposing a model in MOND that suppos- 
edly fits the available data. Here we will examine the 
particular polytropic model proposed by Sll and extend 
the analysis to search for the MOND polytropic model 
that best fits the data on NCG 2419. 

2. THE POLYTROPES 

The model proposed by Sll is a polytrope with poly- 
tropic index n = 10, central velocity dispersion of Co = 
7.5 kms -1 , and central density of po = 35 M Q pc 3 . These 
parameters are related via the polytropic equation: 



«?=cg(p/po) (1/w) 



(1) 



where v 2 is the square of the radial velocity dispersion, 
and p the density. In Paper I, we found that a signif- 
icant orbital anisotropy was necessary to reproduce the 
cluster satisfactorily. The model proposed by Sll also 
include d anisotropy, modeled via th e Osipkov-Merritt re- 
lation ( |Osipkov||1979[ |Merritt||1985 ) 
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The 



where /3 = 1 — Vg/v 2 is the anisotropy parameter. 

variable v 2 is of course the square of the tangential (one- 
dimensional) velocity dispersion. As r a — » oo, the mod- 
els become isotropic. The value adopted by Sll for the 
anisotropy radius was r a — 18 pc. Assuming that the 
cluster is spherical and static, the kinematics and struc- 
ture must obey the spherical Jeans equation: 
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where r is radial distance, and g is the gravitational accel- 
eration. It is straightforward to integrate this system of 
equations numerically, using a simple Euler scheme start- 
ing from r — 0. This procedure works both for Newto- 
nian dynamics and MOND, though to simulate the latter 
the acceleration g needs to be altered according to the 
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Figure 1. Comparison of the surface brightness (top panel) and 
line of sight velocity dispersion (bottom panel) data from Paper I 
to the Sll polytropic model (in MOND). The error bars mark 1- 
a uncertainties. The surface brightness measurements are derived 
from starcounts in annular regions containing between 104 and 910 
stars, whereas the velocity dispersions are calculated from 27 stars 
per bin (except for the outermost bin which has 31 stars). The 
integrated mass of the model is 7.7 X 10 5 Mq, and the mass to light 
ratio has been adjusted (to M/L = 1.8) to obtain the best x 2 fit 
to the surface brightness data. Nevertheless, the fit in panel (a) is 
exceedingly poor, having x 2 = 1123; this can be seen more clearly 
in the middle panel, where we have subtracted the model from 
panel (a). The cross markers show the values of model integrated 
over the same radial annuli as the data, taking into account the 
window function of the imaging survey. Clearly, this particular 
model is completely unacceptable. 



MOND prescription: 



gv(g/ a o) = 9n 



(4) 



where g^ is the corresponding Newtonian acceleration 
that would result from the mass distribution given by 
p(r), while p is the MOND interpolation function. As 
Sll, we take p(x) = x/y/l + x 2 , with cto = 10 -8 cms -2 . 

Figure [l] shows the surface brightness profile and line 
of sight velocity dispersion profile resulting from the pa- 
rameters chosen by Sll. The discrepancy with the ob- 
servations is extremely large. From the 15 points in the 
surface brightness profile alone we obtain a chi-squared 
statistic of x 2 = 1123. While the Sll model is thus 
quantitatively completely excluded, it is very interesting 
nevertheless to examine whether other parameter values 
could provide an acceptable fit. To this end we used a 
general-purpose Markov-Chain Monte Carlo fitting algo- 
rithm to search the space of solutions of the polytropic 
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Figure 2. As Figure [T] but for the most likely stable polytropic 
model fit using the Markov-Chain Monte Carlo algorithm outlined 
in the text. MOND gravity is assumed. As before, the mass to 
light ratio has been adjusted to obtain the best \ 2 fit to the sur- 
face brightness data. We require M/L = 1.9 given the total model 
mass of 7.7 X 10 Mq. However, this model provides a substan- 
tially worse fit to the data than the best Michie model reported in 
Paper I (see Figure [3j: the surface brightness model is a factor of 
81 less likely, and theline of sight dispersion profile is a factor of 62 
less likely. Note that the velocity dispersion data displayed in the 
bottom panel are shown only to guide the eye, as the likelihoods 
are calculated directly from the individual line of sight velocity 
measurements, with their corresponding (individual) uncertainty 
estimates. 

models described above. The input parameters are n, 
Co, po and r a , and the code sweeps through the solutions 
trying to find the most likely model. 

The (log) likelihood is calculated in a similar way to 
Paper I, by adding the logarithm of the likelihood of 
the models given the surface brightness measurements to 
the logarithm of the likelihood of the models given the 
line of sight velocity measurements of individual stars. 
The surface brightness data are derived from starcounts 
measured in annuli over the surveyed area. As detailed in 
Paper I, the HST and Subaru imaging of the cluster has 
significant gaps, so for any accurate model comparison it 
is necessary to take into account the missing areas. The 
window function of the imaging survey was therefore also 
applied to the models, and the models were integrated 
over exactly the same annuli as the data. 

The surface brightness data, by their very nature, are 
binned averages. However, the kinematic data have much 
more discriminating power if they are left unbinned as 
individual line of sight velocity measurements. The like- 



Figure 3. As Figurcp] but showing the best Michie model derived 
in Paper I, assuming Newtonian gravity (#17 from that contribu- 
tion). 

lihood of the model can then be calculated as the prod- 
uct of the likelihood of each data point (thus for these 
likelihood calculations we do not use the binned velocity 
dispersion estimates displayed on the bottom panels of 
Figures Tj(3 which are shown just to guide the eye). For 
the dynamical Michie models discussed in Paper I, we 
integrated the distribution function to derive the line of 
sight velocity distribution as a function of radius. How- 
ever, for the polytropic models analyzed here we did not 
attempt to derive the distribution function, as we judged 
this an unnecessary complication. The reason for this 
is that the observed line of sight velocity distribution 
closely resembles a Gaussian distribution at all radii (see 
Paper I, Figure 14). We therefore assume that the radial 
and tangential velocity distributions are Gaussian, and 
expect this assumption to be a reasonable approxima- 
tion to reality. Furthermore, it is unlikely that any small 
deviations from Gaussianity in the true line of sight ve- 
locity distributions will have any significant bearing on 
the results discussed below. It is straightforward then to 
integrate the polytropic model along the line of sight to 
obtain the predicted velocity dispersion at the projected 
radius of each kinematic datum. The full kinematic data 
set (samples A and B of Paper I) are included in the anal- 
ysis, since this gives the most favorable case for MOND. 
The most likely model found in this manner has 
7.7kms _1 , n = 17.3, r a — 11. Ope and 



Po 



52.7M pc _d . By starting the MCMC algorithm from 
several different initial parameter combinations, we ver- 
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ified that this corresponds to a global likelihood maxi- 
mum. However, this fit did not take into account whether 
the resulting best model is physically plausible. We 
checked that the models obey the Global Density Slope- 
Anisotropy Inequality (Ciotti & Morganti 2010), and 
checked also for stability by calculating ^oi/, a variant 



of the Fridman-Polyache nko-Shukhman parameter (Frid 
man & Poliachenko 19841, defined as twice the ratio of 
radial to tangential kinetic energy within the ha lf-mass 
radius (so th at £hai / = 1 for isotropic systems) . Nipoti 
et al. (20111 have proposed ^haif a s a stability indica- 
tor for spherical stellar systems in MOND, as the maxi- 
mum value for stability £,haif,s is only weakly dependent 
on the density distribution and on the internal accelera- 
tion of the system: in particular, for MOND models of 
NGC 2419 in Paper I we found {half,, ~ 1-4 - 1.5. The 
best model above has £half — 1-53, slightly beyond the 
stable region. Implementing a prior that forces £half < 
1.5 in the fitting procedure, yields a model with parame- 
ters close to the previous best-fit, having cq — 7.9 kms" 



n = 17.6, r a = 11.5 pc and po = 58.0 M Q pc 
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which 



has a total mass of 7.7 x 10° M Q . This best (stable) 
polytropic model is displayed in Figure [2] The residu- 
als with respect to the surface brightness profile are very 
much smaller than for the Sll model, and a first visual 
impression is that this is a reasonable model of the clus- 
ter. Nevertheless, the likelihood of this fit is a factor 
of A = 1/5035 lower than the best Michie model fit in 
Paper I (and shown in Figure [3]): the surface brightness 
and kinematic data contribute a factor of 81 and 62 to 
this likelihood ratio, respectiveljj^] In this comparison 
we have included an 8% component of binaries to the 
Newtonian model, as fit in Paper I; this has the effect of 
producing low-level "wings" to the velocity distribution. 
For the MOND polytropic model, the highest likelihood 
occurs with 0% binaries. 

Both the polytrope and Michie models have the same 
number of parameters: three structural parameters, plus 
an anisotropy radius, plus a fitted mass to light ratio 
and a velocity zero-point. Given that the fits to both 
models have an identical number of degrees of freedom, 
the quant ity —2 In A w ill be asymptotically distributed 
as x 2 (l) (James 2006). This allows us to exclude the 
best MOND polytrope at the 99.996% confidence level. 

While we judge the likelihood approach above to pro- 
vide the strongest and most reliable method for model 
comparison, some readers may prefer the traditional fre- 
quentist \ 2 hypothesis test. In the present context, the 
X 2 method has the disadvantage that it requires the 
kinematic data to be grouped into radial bins to de- 
rive the velocity dispersion. Since the individual veloc- 
ity measurements come from different radial positions 
and have different uncertainties, much information is 
lost in calculating the sub-sample moments. Further- 
more, the (weighted) velocity dispersion estimates listed 
in Paper I have strongly asymmetric non-Gaussian un- 
certainties. The following \ 2 values are therefore pre- 
sented with these caveats. For the best MOND poly- 
trope, we calculate from the 15 surface brightness data 
points Xsb ~ 22.55, whereas the best Newtonian Michie 
model has Xsb = 13-75. From the 6 velocity dispersion 

7 The highest likelihood fit ignoring the stability criterion has 
A = 1/4538. 



measurements (shown in the bottom panels of Figures [TJ 
[3}, we find for the MOND polytrope xlin = 7 - 88 > and for 
the Newtonian Michie model xL„ = 6.76. There are a 
total of 21 data points. The number of degrees of freedom 
for both modelsrlis 15. The probability of an experiment 
with 15 degrees of freedom yielding \ 2 = 30.43 (the value 
for the MOND polytrope fit) or greater by chance is 1%. 
Hence the % 2 test (which for the reasons stated above 
is a weak statistical test for our kinematic data) allows 
us to reject the best polytropic model in MOND at the 
99% confidence level. In contrast, a value of \ 2 = 20.51 
(the case for the best Newtonian Michie model) occurs by 
chance with 15% probability, which is perfectly accept- 
able. Inspection of Figure [2^l shows that the last point 
in the surface brightness profile contributes significantly 
to the model discrepancy. Although there is no reason 
to doubt the validity of this datum, if we ignore it, the 
MOND polytropic model can still be rejected with 97% 
confidence. 

We note finally that it would have been possible to 
perform all the analysis described above by first de- 
projecting the surface brightness distribution to estimate 
the three-dimensional density p(r), as was suggested to 
us by the anonymous referee; the model comparison 
would then have rested purely on the kinematic mea- 
surements. However, the significant uncertainties in the 
surface brightness profile at large radius mean that there 
is no single unequivocal solution to p(r), and any uncer- 
tainty estimates of this function would involve problems 
of correlated noise. Hence our choice to project the mod- 
els into the space of observables should be viewed as a 
statistically simpler option. 

3. CONCLUSIONS 

In the analysis presented above we have derived the 
highest likelihood polytropic model of NGC 2419, fit- 
ting to the observed kinematics and structure of this 
globular cluster assuming MOND and allowing also for 
the possibility of anisotropy in the stellar orbits. The 
model is compared to a previously-fitted model in New- 
tonian gravity. NGC 2419 is probably the best target for 
this anal ysis, since its extreme Ga lactocentric distance 
(87.5 kpc, Di Criscienzo et aL]|2011 ) means that it is rel- 
atively unaffected by the external field due to the Milky 
Way (as confirmed by the N-body experiments presented 
in Paper I), while its large mass and luminosity allow us 
to resolve the radial velocity dispersion profile with a 
useful sample of stars. However, we have found that this 
best-fit polytropic model in MOND is a factor of ~ 5000 
less likely than the best Michie model fit in Newtonian 
gravity, and can be rejected with high confidence. 

Sll suggests that a MOND model deviating slightly 
from the polytropic relation might improve the fit. How- 
ever, we would like to point out that in MOND as in 
Newtonian gravity, given the density profile of a spherical 
model, its velocity dispersion profile is fully determined 
by the anisotropy profile (in other words two different 
distribution functions generating the same density profile 
give different velocity dispersion profiles only if the cor- 
responding anisotropy profiles are different). Thus given 
an observed surface brightness profile, for a choice of M/L 
and P{r), there is only one possible velocity dispersion 

8 Here we use the simpler Newtonian model without binaries. 
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profile for each theory of gravity. In this sense, the exper- 
iment undertaken in Paper I, where we searched through 
a vast number of possible stable anisotropics at different 
M/L, is the definitive test to evaluate the relative like- 
lihood of models. The fact that the polytropic models 
in MOND fare badly in comparison to the Newtonian 
Michie model was anticipated by that experiment. Fur- 
thermore, it must be noted that the best MOND model 
found with the MCMC analysis in Paper I, which was 
shown to perform poorly with respect to the best New- 
tonian Michie model (being less likely by a factor of ~ 
350), turns out to be more likely than the Sll model 
(by a factor of ~ 10 238 ), but also than the best MOND 
polytrope studied here (by a factor of ~ 14). It fol- 
lows that Sanders' (2011) suggestion to use MOND poly- 
tropes, though interesting, did not lead to finding MOND 
models of NGC 2419 better than those already consid- 
ered in Paper I. These results reinforce the conclusions 
presented in Paper I: unless NGC 2419 is significantly 
non-spherical (for which there is no evidence from the 
projected measures we have access to) , or has a radially- 
dependent mass to light ratio (which [Dalessandro et al. 
|2008| argue against, based on the distribution of Blue 
Straggler Stars), or is out of equilibrium, it provides a 
very strong challenge to MOND theory. 

The alternative to this deduction, suggested by Sll, is 
that it is unwise to rely on formal statistical tests given 
that our data and physical analysis may be "plagued by 
systematic effects" . While we have made the upmost ef- 
fort to try to understand and take into account the noise 
properties of the instruments employed in this study (see 
Paper I), it is of course conceivable that random and 
systematic uncertainties bias our results. However, no 
such problems have yet been brought to light, and in 
the absence of such evidence (which we are aware is not 
evidence of absence!), we are unwilling to give up the 
objective tool that is statistical inference. 
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